Phase Field Modeling Of Fracture Research Articles

A phase-field model for the brittle fracture of Euler–Bernoulli beams coupling stretching and bending

Damage gradient models seek to simulate fracture mechanics by modulating the material stiffness. Within this framework, a singular scalar field representing damage is commonly utilized to globally decrease elastic energy. However, when considering structural models like beams and plates, this approach often fails to adequately capture important aspects due to the interaction between stretching and bending contributions. In this work, we propose a model for planar Euler–Bernoulli beams that incorporates the following features: firstly, we utilize two phase-field damage parameters to describe material damaging, specifically addressing the ‘erosion’ occurring above and below the original, undamaged beam; secondly, we assume a simple linear dependence of the axial stress field on the thickness coordinate, along with linear dependence on the axial force and bending moment. By appropriately identifying the constitutive response, our model effectively considers the coupling between stretching and bending induced during through-the-thickness damage, demonstrating good agreement with 2D observations.

Modeling dynamic crack growth in quasicrystals: Unraveling the role of phonon-phason coupling

In this work, we employed the phase-field fracture (PFF) model for dynamic brittle fracture in quasicrystals (QCs). The dynamic PFF formulation is derived using both elastodynamics and elasto-hydrodynamic theory. The temporal discretization is achieved using an implicit generalized-α method, which is unconditionally stable for a particular choice of constants. The model incorporates two distinct crack driving forces: (a) the contribution of elastic strain energy from both phonon and phason fields and (b) the contribution solely from the elastic energy associated with the phonon field. To validate the implementation of QCs, we have conducted benchmarking against recent results in the literature, particularly in the context of quasi-static loading. We addressed various paradigmatic case studies to demonstrate the dynamic crack nucleation and propagation in planar 2D decagonal QCs, with a primary emphasis on understanding the interplay between phonon and phason fields. Phason walls, representing low-energy crack paths resulting from atomic rearrangements, play a significant role in crack propagation by introducing an additional energy component. Notably, increasing the phonon-phason coupling constant expedites both crack initiation and propagation. Furthermore, variations in the phonon-phason coupling parameter result in different crack trajectories, observed in both uniaxial and biaxial loading scenarios. The developed code can be downloaded from https://github.com/Hirshikesh/Quasicrystal.git.

On the energy decomposition in variational phase-field models for brittle fracture under multi-axial stress states

AbstractPhase-field models of brittle fracture are typically endowed with a decomposition of the elastic strain energy density in order to realistically describe fracture under multi-axial stress states. In this contribution, we identify the essential requirements for this decomposition to correctly describe both nucleation and propagation of cracks. Discussing the evolution of the elastic domains in the strain and stress spaces as damage evolves, we highlight the links between the nucleation and propagation conditions and the modulation of the elastic energy with the phase-field variable. In light of the identified requirements, we review some of the existing energy decompositions, showcasing their merits and limitations, and conclude that none of them is able to fulfil all requirements. As a partial remedy to this outcome, we propose a new energy decomposition, denoted as star-convex model, which involves a minimal modification of the volumetric-deviatoric decomposition. Predictions of the star-convex model are compared with those of the existing models with different numerical tests encompassing both nucleation and propagation.

Optimizing integration point density for exponential finite element shape functions for phase-field modeling of fracture in functionally graded materials

In the realm of fracture mechanics and phase-field modeling, the utilization of exponential finite element (EFE) shape functions has shown promise in accurately predicting fracture responses, particularly in scenarios with intricate crack propagation paths. However, the computational complexity associated with EFE shape functions, including higher integration schemes and orientation requirements, prompts the need for optimization. This study delves into the critical aspect of integration point density and aims to determine the ideal balance between accuracy and computational efficiency in the context of EFE shape functions. The research focuses on functionally graded materials, which inherently exhibit spatial variations in solution quantities due to their graded material system. Through a systematic investigation, the study aims to identify the optimal number of integration points required to maximize the computational advantages offered by EFE shape functions while minimizing the associated computational resources. By conducting a comprehensive analysis across various loading scenarios, such as tension and shear, this study intends to establish a quantitative relationship between integration point density and prediction accuracy for fracture responses.

Phase-field modeling of fracture for ferromagnetic materials through Maxwell’s equation

Electro-active materials are classified as electrostrictive and piezoelectric materials. They deform under the action of an external electric field. Piezoelectric material, as a special class of active materials, can produce an internal electric field when subjected to mechanical stress or strain. In return, there is the converse piezoelectric response, which expresses the induction of the mechanical deformation in the material when it is subjected to the application of the electric field. This work presents a variational-based computational modeling approach for failure prediction of ferromagnetic materials. In order to solve this problem, a coupling between magnetostriction and mechanics is modeled, then the fracture mechanism in ferromagnetic materials is investigated. Furthermore, the failure mechanics of ferromagnetic materials under the magnetostrictive effects is studied based on a variational phase-field model of fracture. Phase-field fracture is numerically challenging since the energy functional may admit several local minima, imposing the global irreversibility of the fracture field and dependency of regularization parameters related discretization size. Here, the failure behavior of a magnetoelastic solid body is formulated based on the Helmholtz free energy function, in which the strain tensor, the magnetic induction vector, and the crack phase-field are introduced as state variables. This coupled formulation leads to a continuity equation for the magnetic vector potential through well-known Maxwell’s equations. Hence, the energetic crack driving force is governed by the coupled magneto-mechanical effects under the magneto-static state. Several numerical results substantiate our developments.

A computational meshfree RPIM approach for phase-field modeling of brittle fracture

A computational meshfree RPIM approach for phase-field modeling of brittle fracture

Damage Mechanics Challenge: Predictions based on the phase field fracture model

In this work, we describe our contribution to the Purdue-SANDIA-LLNL Damage Mechanics Challenge. The phase field fracture model is adopted to blindly estimate the failure characteristics of the challenge test, an unconventional three-point bending experiment on an additively manufactured rock resembling a type of gypsum. The model is formulated in a variationally consistent fashion, incorporating a volumetric–deviatoric strain energy decomposition, and the numerical implementation adopts a monolithic unconditionally stable solution scheme. Our focus is on providing an efficient and simple yet rigorous approach capable of delivering accurate predictions based solely on physical parameters. Model inputs are Young’s modulus E, Poisson’s ratio ν, toughness Gc and strength σc (as determined by the choice of phase field length scale ℓ). We show that a single mode I three-point bending test is sufficient to calibrate the model, and that the calibrated model can then reliably predict the force versus displacement responses, crack paths and surface crack morphologies of more intricate three-point bending experiments that are inherently mixed-mode. Importantly, our peak load, crack trajectory and crack surface morphology predictions for the challenge test, submitted before the experimental data was released, show a remarkable agreement with experiments. The characteristics of the challenge, and how changes in these can impact the predictive abilities of phase field fracture models, are also discussed.

A material energy–momentum flux-driven phase field fracture mechanics model

This work introduces a novel phase-field solution for simulating and predicting fractures in elastic solids. In this phase-field fracture model, the crack growth is driven by the material energy–momentum flux, or the configurational force, inspired by a null divergence conservation law and its variational theory, as opposed to employing the commonly-used strain energy degradation-based crack-driving force. By doing so, the crack growth or material configuration change is solely determined by the physics-based Eshelby energy–momentum tensor density or the Rice J-integral flux, which aligns with the energy release criterion of the celebrated Griffith–Eshelby–Rice (GER) theory. The proposed method not only preserves the integrity of the balance of linear momentum principle at the crack tip region but also provides a correct stress asymptotic field in front of the crack tip, while accurately capturing the crack growth.The key idea of this approach is to use the monotonically increasing and irreversible phase field as the marker of the configurational change and form a phase field variational principle that couples the strain energy density with the configurational strain energy density. By doing so, we only damage the displacement field as a form of diffused crack, as opposed to damaging the stress field in the existing phase-field method. Moreover, from numerical simulations, we discovered two different forms of energy release rates: one drives crack growth by enforcing the crack-tip traction loading condition, and one drives crack growth by enforcing the crack-tip deformation condition, demonstrating the ability of the proposed method to naturally capture two distinct modes of fractures due to two different near-crack tip loading conditions. These advancements and contributions reveal previously unknown insights into fracture mechanics.

An efficient and robust staggered scheme based on adaptive time field for phase field fracture model

The standard staggered scheme in phase field model is computationally slow, so that a large amount of iterations are required. In order to improve the computational efficiency, a novel efficient and robust solution scheme based on adaptive time field is proposed in this paper. This scheme is inspired by the adaptive mesh method for space field, and an error indicator based on phase field and history field variables is established, where the time field is coupled with the phase field and energy field in the space field, so that the time field changes adaptively with the computational iterative process, which improves the convergence rate and reduces the time cost. The numerical results show that the presented novel scheme not only significantly improves the computational efficiency, but also has accurate computational results. Compared with the standard staggered scheme, the time cost of the presented scheme is less than 4% of that of the standard scheme. At the same time, the computational efficiency of the presented solver is faster than the current monolithic BFGS scheme, and it has excellent convergence characteristics.

An elastoplastic phase-field model for quasi-static fracture of nickel-based super-alloys

In the present study, an elastoplastic phase-field model of quasi-static fracture in ductile materials is proposed in the variational framework for J2 plasticity with isotropic hardening, which is suitable to describe the quasi-static behavior of metals as investigated in the performed experiments. These contributions include: (1) the free energy functions for coupling elastic response, plastic yielding and damage evolution are established. (2) The new elastic and plastic energy degradation functions are constructed to quantitatively describe the relationship between energy release and phase-field evolution of elastoplastic materials. (3) Damage evolution and plastic yielding criteria are derived. (4) From a numerical point of view, we derive the governing equations and the corresponding weak forms and the overall solution procedure for the phase-field model is given via the use of a return-mapping algorithm. This phase-field model was validated by a series of tensile experiments on Inconel 718 nickel-based super-alloys standard specimens. In order to compare the simulation results with the experimental results more comprehensively, the digital image correlation (DIC) technique is applied to experimentally investigate the specimen deformation information. In addition, to verify the potential of the model to capture complex cracks, we performed Nooru-Mohamed tests. The numerical simulation results are in good agreements with the results of previous experimental work.

The phase-field fracture model enriched by interpolation cover functions for brittle fracture problems

In this study, interpolation cover functions were introduced into the phase-field fracture model (PFM) to improve its convergence and efficiency. The traditional PFM has proven to be successful in simulating brittle fracture problems. However, one disadvantage is that it requires a relatively dense mesh. The PFM model enriched by interpolation cover functions (EPFM) was used for numerical fracture analysis under a relatively coarse mesh. The rank deficiency problem caused by the introduction of interpolation cover functions was resolved using the rank deficiency repair technique. Higher gradients were captured by introducing interpolation cover functions without changing the mesh structure. The numerical examples showed the improvement in convergence. The high efficiency of the proposed EPFM was demonstrated through a series of numerical examples. The numerical results of the EPFM and PFM under different length-scale parameters and meshes were also compared. The results showed that the proposed EPFM is an excellent alternative for evaluating engineering fracture problems.

Rate- and temperature-dependent ductile-to-brittle fracture transition: Experimental investigation and phase-field analysis for toffee

The mechanical behaviour of many materials, including polymers or natural materials, significantly depends on the rate of deformation. As a consequence, a rate-dependent ductile-to-brittle fracture transition may be observed. For toffee-like caramel, this effect is particularly pronounced. At room temperature, this confectionery may be extensively deformed at low strain rates, while it can behave highly brittle when the rate of deformation is raised. Likewise, the material behaviour does significantly depend on temperature, and even a slight cooling may cause a significant embrittlement.In this work, a thorough experimental investigation of the rate-dependent deformation and fracture behaviour is presented. In addition, the influence of temperature on the material response is studied. The experimental results form the basis for a phase-field modelling of fracture. In order to derive the governing equations of the model, an incremental variational principle is introduced. By means of the validated model, an analysis of the experimentally observed ductile-to-brittle fracture transition is performed. In particular, the coupling between the highly dissipative deformation behaviour of the bulk material and the rate-dependent fracture resistance is discussed.

A phase-field fracture model for fatigue behavior in fiber-reinforced composites

A new framework based on the phase field method is developed to model the fatigue fracture behavior of fiber-reinforced composite laminates. The growth path of matrix cracks can be simulated more accurately compared to the classical model. Further, anisotropic fatigue accumulation effects are considered in this model. Comparisons between simulations by the phase field model with uniform fatigue accumulation effects and the proposed method are encouraging. The new framework proposed can more naturally capture different types of crack growth in fiber- reinforced composites. The reliability and potential of this framework to model fatigue behavior for fiber-reinforced composite laminates are demonstrated by investigating several numerical examples. Moreover, the finite element implementation utilizes a simpler and robust approach via the ABAQUS user-defined subroutine UMAT and UMATHT, including the so called AT1, AT2 and PF-CZM models and two solution strategies.

A phase-field fatigue fracture model considering the thickness effect

Modern aero-engine turbine structures often have variable-thickness geometries (e.g., turbine structures with complex cooling systems) that pose a significant challenge to crack analysis. The phase-field fracture method is a numerical method that has recently received considerable attention. However, phase-field fatigue fracture models that consider the effect of thickness are still lacking. To address this problem, this paper proposes a modification of the fracture toughness in the phase-field fatigue model by introducing a relationship function between material fracture toughness and thickness (Gc-B curve). We use reference experimental data to fit the Gc-B curve of GH4169 and predict the fracture toughness of specimens with other thicknesses. The predictions of the proposed model agree well with the experimental results reported in the literature. This achievement contributes significantly to the structural safety design and life prediction of aero-engine components.

Modeling brittle crack propagation for varying critical load levels: a dynamic phase-field approach

AbstractBrittle materials are known for their violent and unpredictable cracking behavior. A behavior which is dictated by a combination of microscopical material defects and the competition between the potential energy of the system and the surface energy of the material. In this study, we present the implementation of a dynamic fracture phase-field model with a new crack driving force into a commercial finite element (FE) solver and examine its behavior using three different tension-compression splits. After validating the implementation, we use the model to investigate its predictive capacity on quasi-statically loaded L-shaped soda-lime glass specimens with varying critical load levels. The dynamic fracture phase-field model predicted similar crack propagation to what was found in the literature for quasi-static and dynamic validation cases. By varying the critical load level for the L-shaped soda-lime glass specimens using the new crack driving force, the model predicted a positive correlation between the initial crack propagation speed and the critical load level, similar to what was seen in the experiments. However, the predicted crack propagation speed decreased quicker than the experimental crack propagation speed. The tension-compression splits had an impact on the predicted crack propagation paths. Overall, the proposed crack driving force used in the dynamic fracture phase-field model seems to capture the relation between critical load and initial crack propagation speed and thus enables crack predictions for specimens of varying strength.

Fracture behaviors of nacre-like composites via phase-field fracture modeling

Nacre-inspired composites, composed of hard and soft phases, exhibit exceptional combinations of stiffness, strength, and fracture toughness. However, the fracture behaviors of these composites are significantly influenced by the geometric factors and material properties of their constituents. In this study, we systematically investigate the fracture behaviors of nacre-like composites by employing a phase-field fracture model considering the effects of oblique angle between loading and platelet directions, a critical factor often overlooked in previous studies, coupled with the effect of other geometric and material parameters, such as the aspect ratio, Young's modulus and fracture toughness ratio between hard and soft phase. Our findings highlight the significant influence of the oblique angle on the mechanical properties and fracture behaviors of these composites. Particularly, it is found the maximum Young’s modulus and strength are obtained when the direction of hard platelet is aligned with the tension direction, while the maximum toughness is achieved at a moderate oblique angle, especially in cases of large aspect ratios, small Young's modulus and fracture toughness ratios between the hard and soft phases. Besides, the maximum fracture toughness is more likely to be obtained at a small aspect ratio and large Young's modulus ratio between hard and soft phases. The insights gained from our study would offer valuable guidance for the design and practical application of nacre-mimetic composites.

Strengthening mechanism of steel fiber in UHPC: A new fracture phase field model

Strengthening mechanism of steel fiber in UHPC: A new fracture phase field model

Phase-field modeling of anisotropic crack propagation based on higher-order nonlocal operator theory

This paper presents a novel higher-order nonlocal operator theory for the phase-field modeling of brittle fracture in anisotropic materials. Incorporating higher order nonlocal operators can enhance the accuracy of the phase-field model by effectively capturing long-range interactions that hold significance in numerous materials. The reproducing kernel particle method is employed to derive a nonlocal differential operator to enhance computational stability and accuracy. Moreover, the proposed method eliminates the need for direct computation of derivatives of the modified kernel function, which avoids the calculation of moment matrix derivatives and improves computational efficiency. The phase-field modeling of polycrystalline materials, considering the anisotropic fracture resistance of each grain, is implemented using this numerical framework. The present method is able to capture different scenarios intergranular and transgranular crack propagation patterns in polycrystalline materials. The proposed method involves a detailed representation of the complex process of crack initiation and propagation in 2D and 3D models of polycrystalline materials.

Adaptive meshfree method for fourth-order phase-field model of fracture using consistent integration schemes

An adaptive fourth-order phase-field model is proposed using consistent element-free Galerkin method (EFG). Convenient construction of high-order interpolation functions in EFG method is fully utilized. Requirement of C1 continuity of shape functions is satisfied in fourth-order phase-field model. Two different consistent integration schemes with high efficiency and high accuracy, namely quadratically consistent one-point integration scheme and three-point integration scheme, are employed to evaluate the fourth-order phase-field equation and displacement equation, respectively. To reduce computational cost, the adaptive strategy of local background mesh refinement is established on basis of strain energy history and phase field. Background mesh is refined via the insertion of nodes at the midpoints of each side of the integration element. Comparing with EFG method, computational accuracy and efficiency are enhanced by the proposed method. Load-displacement response, the smoothness of the resulting stress and crack paths are predicted well.

Evaluating tensile strength of cement paste using multiscale modeling and in-situ splitting tests with micro-CT

This study investigated a strategy to simulate splitting tests to evaluate the tensile strength of cement paste within a multiscale framework, and validated the approach with in-situ micro-CT measurements. Phase-field fracture model was used to simulate the splitting test using virtual microstructures reconstructed from the macroscale micro-CT measurements. Among the modeling parameters, the range of the diffusive crack width, which is an internal parameter for phase-field fracture model, was further investigated. The cracks propagated through the weaker solid phase, and the predicted strengths and standard deviations were comparable to those obtained from the experiments. It was found that the diffusive crack width should be smaller than the characteristic length of the weaker phase to predict the crack pattern accurately. These findings confirm the effectiveness of the proposed multiscale modeling approach in predicting the tensile strength of cement paste.